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We study the finite-temperature behavior of the condensate in the Landau gauge of SU{2) 
Yang-Mills theory on the lattice in a wide range of temperatures. The asymmetry between the elec- 
tric (temporal) and magnetic (spatial) components of this unconventional dimension-2 condensate 
is a convenient ultraviolet-finite quantity which possesses, as we demonstrate, unexpected proper- 
ties. The low-temperature behavior of the condensate asymmetry suggests that the mass of the 
lowest thermal excitation in the condensate is unexpectedly low, about 200 MeV , which is much 
smaller than the glueball mass. The asymmetry is peaking at the phase transition, becoming a 
monotonically decreasing function in the deconfinement phase. A symmetric point is reached in 
the deconfinement phase at a temperature approximately equal twice the critical temperature. The 
behavior of the electric-magnetic asymmetry of the condensate separates the phase diagram of Yang- 
Mills theory into three regions. We suggest that these regions are associated with the condensed, 
liquid and gaseous states of the confining gluonic objects, the Abelian monopoles. 

PACS numbers: 12.38.Aw, 11.15.Ha, 25.75.Nq 
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I. INTRODUCTION 

In a specific theoretical frame, non-perturbative fea- 
tures of QCD are reflected by the existence of vari- 
ous non- vanishing local condensates. The most famous 
vacuum condensates are the dimension-4 gluon conden- 
sate, (0^(0°^)^), and the dimension-3 quark condensate, 
{tl^ij}) ■ The former characterizes the non-perturbative dy- 
namics of strongly interacting gluon fields while the latter 
- in the formal limit of vanishing quark masses - is an 
order parameter corresponding to the spontaneous break- 
ing of chiral symmetry. The local condensates enter the 
QCD sum rules as non-perturbative power corrections 
having great significance for QCD phenomenology 

The dimension-2 condensates represent somewhat un- 
conventional vacuum condensates. Indeed, it is impos- 
sible to construct dimension-2 operator in QCD in a 
gauge-invariant and local manner, simultaneously. One 
has to abandon either the condition of locality or the 
requirement of gauge-invariance. Nonlocal but gauge in- 
variant operators are useless from the point of view of 
the operator product expansion: one cannot relate a non- 
local operator to a Wilson coefficient corresponding to a 
power correction using the dimensionality counting rule. 
On the other hand, the local but gauge-dependent op- 
erators seem to be useless because, as one might think 
naively, these operators cannot contribute to physical 
observables. However, this conclusion seems to be not 
compelling as it was pointed out in Ref. [1, Q . 

The siniplest dimension-2 operator in SU{Nc) gauge 
theory is 
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(1) 
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This operator is not gauge invariant. Therefore its expec- 
tation value can be understood in a two-fold way: one can 
either average the operator over all possible gauge trans- 
formations or one can evaluate the operator at a partic- 
ular point of a gauge orbit. The former choice leads to 
the unwanted non-locality while the latter is a suitable 
option in particular gauges. The Landau gauge is defined 
as the result of minimization of the bulk averaged A^- 
operator with respect to gauge transformations. Thus, 
the extremum of the dimension-2 operator gets a special 
meaning since A^^x) is a local operator in the Landau 
gauge. 

The dimension-2 condensate enters the ultraviolet 
asymptotics of the gluon propagator in the Landau 
gauge as a non-perturbative power- like correction. 
Moreover, this condensate emerges also in QCD in the 
quark propagator and in various vertices in the Landau 
gauge [5j- It is worthwhile to remember that for a long 
time it was assumed that due to asymptotic freedom the 
whole physics at short distances must be described by 
exclusively perturbative physics. The appearance of the 
dimension-2 operator in the ultraviolet regime is a sig- 
nature of mixing of non-perturbative and perturbative 
features of QCD: non-perturbative effects are emerging 
at very short distances. 

Besides the purely theoretical issue of mixing ul- 
traviolet and infrared physics, the appearance of the 
dimension-2 condensate plays a distinguished role in 
QCD phenomenology because it is associated with non- 
standard power corrections 0, d, Q . This motivates the 
wide interest in the subject, such that the A^-condensate 
was intensively studied both numerically and analyti- 
cally. Numerical simulations of lattice SU (3) Yang-Mills 
theory at zero temperature indicate that the conden- 
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sate is, in fact, a large quantity 

= [1.64(15) GeVf. (2) 

The energy scale of the dimension-2 condensate ([2|) is 
of the order of the glueball mass [1, 0| 

TOO++ = 3.52(12) v/cr(r = 0) = 1.55(5) GeV, (3) 

if one takes the phenomenologically accepted value a{T — 
0) — [440 MeV]^ for the string tension at zero temper- 
ature. The coincidence does not look accidental. Of 
course, the appearance of the dimension-2 condensate 
and the mass gap generation in the non-Abelian gauge 
theory are both of non-perturbative nature. 

Another non-perturbative phenomenon in non-Abelian 
gauge theory is the confinement of color which is often 
linked with the mass gap generation but the relevant scale 
is Aqcd ~ 200 MeV. A relation between the dimension- 
2 condensate and color confinement may also be existent 
although the very reason is not clear at present. A good 
playground to study the interrelations between the dif- 
ferent mentioned non-perturbative features would be a 
suitably chosen effective (toy) model. The confinement 
of electric charges and the mass gap generation is un- 
derstood, for example, in the framework of an Abelian 
gauge model with a compact gauge field, which can be 
considered as a certain limit of the 5*0(3) Yang-Mills- 
Higgs model (often called Georgi-Glashow model). The 
compactness of the gauge field is related to presence of 
non-perturbative objects, magnetic monopoles, the dy- 
namics of which leads both to confinement and to mass 
gap generation [lo| . 

The four-dimensional version of the compact Abelian 
gauge model (sometimes called "compact electrodynam- 
ics" or cU{l) gauge theory) contains two phases - the 
confinement phase and the deconfinement phase - sep- 
arated by a phase transition. At zero temperature this 
model is especially interesting since the phases of the 
model are characterized by a single parameter, the gauge 
coupling g. The confinement phase is located at strong 
coupling, g > ~ 1, while the deconfinement regime 
is associated with weak coupling. The confinement phe- 
nomenon comes along with the mass gap generation while 
in the deconfinement phase the mass gap shrinks to zero. 
Both mass gap and charge confinement have the same, 
well understood [ll[ roots in the condensation of the 
Abelian monopoles (for a review see Refs. [H, [3|). The 
natural question to ask is "does an A^-condensate appear 
in the compact Abelian model?" If the answer is posi- 
tive, what can we say about a possible relation between 
confinement/mass gap generation and the emergence of 
the condensate? These questions were addressed in 
Ref. 0, where it was found that the non-perturbative 
part of the condensate in the Landau gauge is a very 
good order parameter for the confinement-deconfinement 
phase transition: this dimension-2 condensate is non- 
vanishing in the confinement phase and vanishes in the 
deconfinement phase. In the compact Abelian model the 



link between all three phenomena - confinement, mass 
gap generation and emergence of the dimension-2 con- 
densate is obvious: the primary reason of all these phe- 
nomena is monopole condensation. 

Coming back to the Yang-Mills theory, one finds an 
essential difference between this theory and the just dis- 
cussed Abelian model: the transition in the pure Yang- 
Mills theory is driven by thermal fiuctuations and, there- 
fore, it happens at finite temperature. On the contrary, 
in the mentioned simulations of the compact Abelian 
gauge theory the transition is purely quantum, and ther- 
mal fluctuations are not involved. In finite tempera- 
ture Yang-Mills theory (and in QCD as well), instead 
of having one dimension-2 condensate ([T]) one needs to 
define two types of it, corresponding to time- like ("elec- 
tric") and to space-like ("magnetic") gauge bosons, sep- 
arately. In the four-dimensional Euclidean space-time 
corresponding to the imaginary-time formulation of the 
theory at T ^ 0, one defines the electric and magnetic 
contributions, respectively: 

Al = l-TTAi{x)A4x), 

< = ^Tr^A,(x)A,(:r), (4) 

^ i—l 

such that the full dimension-2 condensate at nonzero 
temperature is the sum of both, 

(g^A^) = (g^Al) + (g^Al,) . (5) 

In Eq. ^ we have not yet divided the magnetic contri- 
bution by the number of spatial dimensions. This would 
make {g'^A\) and {g'^A\j) the dimension-2 condensate 
per one Lorentz component (or per three Lorentz com- 
ponents) of the electric (magnetic) gluons in the Landau 
gauge. We refine the conclusion of Ref. [l3| that the 
A^ condensate observed at zero temperature is consis- 
tent with a vanishing condensate in the deconfinement 
phase. In fact, the condensate at finite temperature is 
characterized by the electric and magnetic components 
which, as we show, are quite nontrivial. To our knowl- 
edge, the eventual difference between the two simplest 
dimension-2 condensates was not considered so far. In 
this paper we fill this gap. 

Before proceeding further one should comment about 
the ultraviolet divergences of the dimension-2 conden- 
sate. In compact electrodynamics [3] at vanishing tem- 
perature the condensate contains perturbative and non- 
perturbative parts, respectively, 

{g'A^)^{g^A^)^^^^, + {g^A^)^^. (6) 

The perturbative part is quadratically divergent whereas 
the non-perturbative part is finite. The first one can be 
calculated trivially in non-compact QED, the free the- 
ory of photons. It is the non-perturbative part which 
serves as the order parameter of the non-thermal phase 
transition in compact electrodynamics. 
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In QCD the situation is similar but not so sim- 
ple 0, [isl . [iGj . A linear decomposition (O works as 
well up to renormalization-related logarithmic correc- 
tions. Then in Eq. ® {g'^A^)^^^^^ oc A^v and {g'^A'^)j^p oc 
^QCD' where Aqcd defines a typical QCD energy scale - 
supplemented with a large prefactor according to Eq. 
- while Auv oc 1/a is an ultraviolet cutoff, which in the 
case of lattice calculations is inversely proportional to the 
lattice spacing a. The calculation of the non-perturbative 
part of the condensate requires a renormalization which 
was implemented, for example, in Ref. 7]. 

At finite temperature the decomposition © should 
hold as well. It is well known that the finite-temperature 
corrections to physical observables do not contain ultra- 
violet divergences. However, the temperature corrections 
can be both of perturbative and non-perturbative nature 
such that they can affect both the perturbative and non- 
perturbative parts of the decomposition 

In our numerical analysis we do not discriminate 
between perturbative and non-perturbative parts. Of 
course, we distinguish between corrections received by 
the spatial (magnetic) and temporal (electric) gluons. On 
general grounds we write the decomposition per Lorentz 
component in the form 
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{g'A\ + {g'Alh 



(7) 



asymmetry ([9]) is a finite-valued quantity in the ultravi- 
olet regime. 

We have calculated the asymmetry ([9]) using numeri- 
cal simulations of Yang-Mills lattice gauge theory in the 
Landau gauge. Obviously, at zero temperature (symmet- 
ric lattice) all Lorentz components of the gauge field Afj_ 
contribute equally to the A"^ condensate because of the 
approximate 0(4) rotational symmetry (actually H{A) 
hypercubic symmetry) satisfied by the Euclidean lattice. 
At finite temperature this rotational symmetry is broken 
down to (approximate) 0(3) spatial symmetry. Space 
and time coordinates are no more equivalent since, in 
the imaginary time formalism used to study the field in a 
thermodynamic equilibrium state, a finite temperature T 
is implemented via compactification of the (imaginary) 
time direction to a circle of length Lt = 1/T, while the 
other (spatial) directions Lg are still infinite, or at least 
Lg ^ Lt- Thus, the fluctuations of the temporal (elec- 
tric) A4 and spatial (magnetic) Ai {i — 1,2,3) compo- 
nents of the gluon fields must be different in general. 

The structure of this article is as follows. In the fol- 
lowing Sect.|TT]we describe the theoretical expectations. 
In Sect, mil we report our lattice simulations for SU{2) 
pure gauge theory. Sect. IIVI contains a discussion relat- 
ing our findings to monopoles and confinement. We draw 
conclusions and define further routes in Sect. fVl 



{g'A')M = l{g'A^), + {g'Al,)T, (8) 

where {g'^A^)^ is the dimension-2 condensate at zero 
temperature while finite-temperature 
correction to the electric (magnetic) dimension-2 con- 
densate. These thermal corrections are ultraviolet-finite 
at finite T and are vanishing at T = 0. It is the zero- 
temperature condensate {g'^A^)^ that contains a piece 
which is quadratically divergent in the ultraviolet. In 
the limit of zero temperature magnetic and electric com- 
ponents are equal and expressed naturally via the zero- 
temperature condensate. In this limit the sum rule ([5]) is 
also naturally restored. 

We suggest to concentrate on the thermal corrections 
to the condensates, {g'^A\)T and {g'^A\j-)T, as probes of 
the thermal activity of the electric and magnetic compo- 
nents, respectively, of the gluonic medium. In this pa- 
per we compute the electric-magnetic asymmetry of the 
dimension-2 condensate 



A^2(T) 



{g'A^E) 

.2 /i2 
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{9'^%)T~^^{g'Ai,)r 



{g'Ai,) 
1 

3' 



(g'Ai)^ \ y: (g'A^)^ . 



(9) 



The asymmetry (|9]) plays a special role since the quadrat- 
ically divergent zero-temperature components of both 
electric and magnetic condensates cancel out. Thus the 



II. THEORETICAL EXPECTATIONS 

The low- and high-temperature asymptotics of the 
electric-magnetic asymmetry of the dimension-2 conden- 
sate A^2 can be guessed from general arguments: 



Aa2{T) cx 



e 

y2 



T:$>Tc. 



(10) 



It is tempting to equate the mass parameter rrigi with the 
mass of the lowest glueball: 



[naive expectation] . (11) 



In Eq. pop a polynomial prefactor at low temperatures 
and possible logarithmic corrections at high tempera- 
tures are omitted. The high temperature asymptotics 
in Eq. (jlOp is determined for dimensional arguments. In 
the result of lattice simulations we will show that the 
dimensional argument is correct, as expected. 

Common wisdom says that the low-temperature 
asymptotics of any thermodynamic system is determined 
by properties of the lowest excitation. The lowest excita- 
tion in the SU{Nc) gauge theory is a glueball, which has 
a non-zero mass due to the phenomenon of the mass gap 
generation. According to Ref. [81], in the SU{3) gauge 
theory 



mo++ = 3.52(12)v/cr(T = 0) . 
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According to the Bielefeld group, Ref. [13, llS] the finite 
temperature phase transition in pure SU{S) gluodynam- 
ics happens at 



such that 



Tc = 0.629{3)^/a{T = 0) 



mo++ ~ 5.6 Tc 



Given that mQ++ ^ Tc the low-temperature asymptotics 
should hold true not only at T Tc, but also at temper- 
atures close to the phase transition temperature. As we 
will see below, in Yang-Mills theory the exponential form 
of the low-temperature asymptotics is correct while, un- 
expectedly, the "natural" identification of mass scale (fTT|) 
turns out to be wrong. Moreover, the low-temperature 
asymptotics (ITU)) of the asymmetry turns out to be also 
incorrect in certain simple models indicating that the ex- 
ponential low-temperature asymptotics (jlOp is not valid 
in a general case. 

We begin the discussion of the asymmetry in the sim- 
plest case of free photodynamics and then for the Abelian 
Higgs model. After that we turn to the case of Yang- 
Mills theory. 



A. Example of photodynamics 

Photodynamics is the theory of the free Abelian gauge 
field. Since the theory does not contain any dimensional 
scale and is not able to generate a scale by itself like QCD, 
the electric-magnetic asymmetry of the condensate 
- which is the dimension-2 quantity A^2(T) - must be 
proportional to for dimensional reasons. 

Since the Lagrangian of photodynamics is quadratic, 



•^phot — —F, 



F^, = d,A. - d,A^ , (12) 



it is easy to calculate the photon correlation function 
(propagator) in momentum space 



(13) 



where p = (p,P4) is the 4-momentum, while the relation 
between the photon fields in coordinate and momentum 
spaces is given by Fourier transformation. At zero tem- 
perature the relation is 



d^p 



e^iP^-)A^ip), 



(14) 



while at finite temperature it is given by the formula 

'^'^ TVe'(P'-^+*""-M^(p,c.„), (15) 



(16) 



where 



uJn = 2'KnT 



are the Matsubara frequencies. 

Due to absence of interactions the form of (zero tem- 
perature) photon propagator 



(17) 



D,,{p)^D{p'){6,,-^-^ 



is also valid at finite temperature. The propagator (jl7p 
is parameterized by a single propagator function (here 

p^ = 



p. 



(18) 



The asymmetry of the dimension-2 condensate in pho- 
todynamics is calculated in Appendix \X[ 



y^free 



(T) = 



OO 

^/d..V.(p){^-i[i + M.)]} 
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(19) 



As the definition of the asymmetry we took Eq. ^ 
with one-component gauge field, Ng — 1. We have also 
omitted the overall coefficient, electric charge squared, 
— e'^, in the definition Q in order to simplify the ex- 
pressions below. This coefficient can easily be restored: it 
should enter all the analytic results of this Section for the 
asymmetry as just the proportionality coefficient. The 
proportionality of the asymmetry to the squared temper- 
ature in Eq. (|19p is quite obvious due to dimensionality 
reasons. The fact that the asymmetry ^T9\\ is negative 
tells us that in the absence of the interactions the fluctu- 
ations of the magnetic (spatial) photons are dominating 
the fluctuations of the electric (temporal) photons. 



B. Example of Abelian Higgs model 

Now we turn to a more complicated case, adding a 
charged scalar field to photodynamics. This system is de- 
scribed by the Abelian Higgs model with the Lagrangian 



HM - \f',. + \id^ + ^eA^)^\^ + A(|0p - rj^ , (20) 



and controlled by the gauge coupling e and the quartic 
coupling A. The overall dimensional scale is fixed by the 
position rj of the minimum of the potential for the Higgs 
field 0. 

In the Higgs phase the photon has a mass m — erj due 
to spontaneous symmetry breaking caused by the con- 
densation of the scalar Higgs field. At zero temperature 
the photon propagator is described by Eq. (|17p with the 
single propagator function 



1 



p^ 



(21) 



In the case of the London limit, A — > oo, with the Higgs 
field being infinitely massive, the latter is not thermally 
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excited if the system is subject to finite temperatures. As 
the result, the Higgs loops do not contribute to the pho- 
ton polarization tensor. The difference between spatially 
longitudinal and spatially transverse photons is absent in 
this case as they share the same propagator function pT|) . 
Then the Lorentz structure at finite temperature of the 
gauge propagator stays the same as in Eq. (|17p. 

The asymmetry for the London limit is calculated in 
Appendix [B] 



e(T,TO)-£(T,TO = 0) 



-il](T,m) , 



(22) 



where 



s{T,m) = J ^Vp2+m2 j^(^p2+^2) (23) 



and 

s(r, m) 



(2^ vFT^ 



MVp' + m^) ■ (24) 



In Figure [T] we show the asymmetry (pS)) divided by 
the temperature squared, T^, as a function of the tem- 
perature T expressed in units of the gauge boson mass m. 




T/m 



FIG. 1: (Color online) The normalized electric-magnetic 
asymmetry of the condensate, A^2(T, m)/r^, as a func- 
tion of the normalized temperature, T/m, in the Abelian 
Higgs model. The horizonal dashed line in the plot shows 
the high-temperature limit (|27p . for T ^ m, which recov- 
ers the case of photodynamics. The inset illustrates the low- 
temperature behavior of the condensate asymmetry 



In the massless limit, 
be calculated exactly: 



0, the asymmetry can 



lim AAHM(r^^)^_Z_, 



(25) 



The comparison of Eq. ([25)) with Eq. (|T9|) shows that 
the -asymmetry is an analytic function of the photon 
mass at m = 0: 



lim A^?^(T,m) 



= A 



AHM 
A2 



{T,0) = Af°'°{T). (26) 



This result is, in fact, not guaranteed from the beginning 
- as we describe this in the Appendices - since the ex- 
actly massless case involves the calculation of a residue 
at the double pole while the massive theory has always 
two isolated single poles in the asymmetry. 

Equation ([25]) corresponds, in fact, to the leading term 
in the asymmetry in the high temperature limit, T ^ 
TO. Supplementing this result with the low-temperature 
expansion we get the leading terms in high- and low- 
temperature limits, respectively: 



J.2 

A^2 (T,m) = 



Aji2 (T, m) 



27r2 
45 TO^ 



(T»to), (27) 

(r<TO), (28) 



where the ellipsis denote subleading contributions. 

As we have expected the difference between the elec- 
tric and magnetic components of the condensate in the 
Abelian Higgs model (AHM) is vanishing in the zero tem- 
perature limit according to Eq. ([28]) . What comes un- 
expected is that at low temperatures the asymmetry is 
suppressed polynomially (by the fourth power of the tem- 
perature, r^) and not exponentially as one would have 
expected from general arguments presented in the begin- 
ning of this Section. Indeed, the Abelian Higgs model is 
a theory with a finite mass gap, and therefore one could 
expect that thermodynamical contributions to any quan- 
tity would be suppressed exponentially in the low tem- 
perature limit by a factor exp{— to/T}. Our analytical 
calculation shows this is not the case. 

Technically, the unusual polynomial behavior 
A^2 (T, to) (X of the asymmetry at low temper- 
atures is due to the massless, m — 0, term in Eq. (|22p . 
This term corresponds to a massless longitudinal degree 
of freedom while all other terms are massive such that 
they are suppressed exponentially as T -C to. The ap- 
pearance of this term can be traced in the integrand (|B2[) 
of the integral representation of the asymmetry (|AII[) 
as this integrand contains a massless pole. The massless 
pole, in turn, appears due to the fact that the propagator 
of the gauge boson p7|) contains a I/p^ term which 
does not give any contribution in explicitly transverse 
gauge-invariant expressions such as, for example, the 
correlator of two field strength tensors. However, the 
A^-propagator does incorporate an infrared l/p^ term, 
which gives rise to a polynomial behavior unless the 
propagator function D{p) is vanishing in the infrared. 
Both in photodynamics and in the Abelian Higgs model 
the photon propagator is either divergent or finite in the 
infrared limit such that the zero temperature limit of 
the asymmetry is not exponentially suppressed contrary 
to our naive expectation ([TI 
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C. Distinguishing longitudinal and transverse 
photons 



III. NUMERICAL RESULTS IN SU(2) 
YANG-MILLS THEORY 



Contrary to the considered examples, at finite temper- 
ature the gluon propagator in the Landau gauge is in gen- 
eral parameterized by two propagator functions. These 
are the transverse (or, "magnetic") propagator Dt and 
the longitudinal (or, "electric") propagator D^. In mo- 
mentum space, 

i?^t(p,P4) = S'''[p;[, Drip^pi) + P^, Z?l(p,P4)] . (29) 

Here and P^^ are, respectively, projectors onto spa - 
tially transverse and spatially longitudinal directions [1^ , 
with 



(1 - 5^i) (1 - S,i) 1 



satisfying the relation 

pT , pL _ p p_x _ 



p2 



(30) 



(31) 



where Pp,y is the standard 0(4)-symmetric projector cor- 
responding to the zero temperature case with 

At finite temperature the propagators Dt and 
are, in general, different from each other. The differ- 
ence between spatially longitudinal and spatially trans- 
verse functions arises due to interactions among the fields 
while the free gauge theory has these functions equal: 
P)^°° = D^l'"^. The interactions which lead to the differ- 
ence between the propagators Dt and Dt may be pertur- 
bative, as, for example, in quantum electrodynamics p^ . 
or these interactions can be of purely nonperturbative 
origin, as, for example in an Abelian gauge theory with 
contains only a compact gauge field [20|. The compact 
Abelian gauge theory possesses nonperturbative topolog- 
ical defects, monopoles, which affect drastically the prop- 
agator properties in the confining phase of the theory. 

In SU{2) Yang-Mills theory the propagators Dt and 
Dt were investigated using both numerical simulations 
on the lattice zl'l and analytical calculations in the con- 
tinuum ^22, 2^, ,24 , 25]. In the limit of vanishingiy low 
temperatures the thermodynamics of Yang-Mills theory 
imposes a certain constraint on the infrared critical ex- 
ponents which characterize the infrared behavior of the 
correlators [1^. 

From Eq. (|29p one can derive the electric-magnetic 
asymmetry of the condensate in terms of the two 
propagators: 



1. 



24 7r3 



-T 



r3p2 



-DT{p,L^n) 



-2D' 



(32) 



This expression should be ultraviolet finite up to a loga- 
rithmic renormalization. 



In the following we report measurements of the 
electric-magnetic asymmetry of the condensate that 
were performed in SU (2) gluodynamics simulations. The 
configurations have been created by means of a heatbath 
Monte Carlo code. The minimal Landau gauge has been 
enforced for every 30-th configuration before measuring 
the . For the gauge fixing we have used the Simulated 
Annealing algorithm 27]. An interval of the gauge tem- 
perature ranging from Tmax = 1-0 to T,nin = 1.0 • lO^'^ 
has been traversed within 3000 sweeps with linearly de- 
creasing gauge temperature. This "gauge cooling" was 
followed by obligatory over relaxation until the required 
transversality was reached. The stopping criterion was 
max^maxa |(5^A°)(a;)| < 10"^. 

We have evaluated the electric-magnetic asymmetry 
on lattices 16^ x 4, 24^ x 6 and 32^ x 8. In the interval 
[3 e [2.20, 2.95] we have selected a grid containing 51, 36 
and 24 /3-values, respectively. In this way, an interval of 
physical temperatures T S [0.4rc, S.lTc] is covered. In 
the restricted range T G [0.4rc, 2.5rc], systems at nearly 
equal temperatures are realized by lattices with different 
lattice spacings. Close to the deconfining transition, the 
31? volumes are equal to each other with an aspect ratio 
of A^s : -^t = 4 : 1. The simultaneous evaluation of 
the asymmetry provides us with a valuable assessment of 
potentially dangerous finite cutoff effects. This does not 
seem to be a problem at all. As we will see later, data 
from different lattice sizes are smooth and can be fitted 
simultaneously as functions of the physical temperature. 

The number of configurations at each combination of 
lattice size and /?- value was adjusted in such a way to give 
reasonable statistical error bars (typically, the errors are 
of the order of a few percent near the phase transition 
temperature). Close to the transition we needed about 
1500 configurations per /3 value at the smallest lattice, 
and we found about 50 configuration per value of [3 suf- 
ficient at the largest lattice far from the transition. 

All our production measurements have been done with 
one Gribov copy only. We stress that the Simulated An- 
nealing algorithm already shifts the outcome of a sin- 
gle gauge fixing closer to the (unknown) absolute maxi- 
mum of the gauge functional than several repetitions of 
pure overrelaxation could do. We have actually checked 
the Gribov copy dependence at our middle-sized lat- 
tice, 24"^ X 6 for a representative set of four tempera- 
tures corresponding to the confinement region [(3 = 2.30, 
T K, 0.65 Tc), to the deconfinement region (/? = 2.7, 
T K, 2.4 Tc) and to two temperatures close to the phase 
transition, on the confinement (/? = 2.40, T « 0.9 Tc) 
and on the deconfinement (/3 = 2.50, T « 1.25 Tc) sides. 
We took the first ten Gribov copies, Nq = 1 ... 10, by re- 
peating the Simulated Annealing gauge fixing. To form 
the ensemble collecting the best gauge copies for each 
Monte Carlo configuration after Nq repetitions of Sim- 
ulated Annealing we sampled the "currently best" copy 
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among the preceding Nq gauge-fixing trials. We found 
that - within our statistical errors ~ the plateau value (for 
the ensembles with Nq — > oo) of the electric-magnetic 
asymmetry is already reached in the ensemble of gauge- 
fixed copies corresponding to Nq — 4 trials. In the de- 
confinement phase the copy dependence is negligible since 
the systematical uncertainty due to the Gribov copy de- 
pendence is much smaller than the statistical errors of our 
calculations. In the confinement region the uncertainty 
- calculated as the relative deviation of the measured 
value of the electric-magnetic asymmetry after Nq trials 
from the Nq — > oo plateau value - is about 2% while in 
the close neighborhood of the phase transition the sys- 
tematic correction to the electric-magnetic asymmetry 
due to Gribov copy dependence may reach 10%. The 
asymmetry is slightly rising with the number of gauge 
copies being under inspection. However, all characteris- 
tic features of the asymmetry discussed in this article are 
unaffected by this Gribov copy dependence. 

The vector potential is extracted from the links as 

gA;{x + fi/2)= Tr[r, (f/,,^ - [/] J /(2 * a)] . (33) 

We express all dimensional quantities in units of the 
critical temperature. According to Ref. the phase 
transition - which is of the second order in pure SU (2) 
gauge theory - happens at the critical temperature 



Tc = 0.694(18) V'7(r = 0) = 305(8) MeV. 



(34) 



The electric-magnetic asymmetry of the A^- 
condensate ^ is shown in Figure [51 One can observe a 



<v, 

< -2 
-4 
-6 



r=0.694(18) a =305(8) MeV 
symmetric point, 7^=2.21(5) J 
/ 




16x4 
24'x6 
32'x8 

■ fit (777=1.5. ..6.2) 



T/T 



FIG. 2: (Color online) The electric-magnetic asymmetry of 
the j4^-condensate ((9]) for SU{2) gauge theory normalized by 
the critical temperature squared as function of T/Tc. The 
high-temperature fit (|39|) with the best fit parameters (|40p 
is shown by the solid line. The horizontal dashed line cor- 
responds to A^2 = while the vertical dashed line marks 
the critical temperature. The symmetric point is explicitly 
indicated. 

very good scaling: the results obtained at various lattice 
sizes describe the same curve when plotted in physical 



units. This fact indicates that the lattice ultraviolet 
artifacts are negligible in our numerical setup. 

Before going into details we notice immediately two 
distinct features of the temperature dependence of the 
condensate asymmetry: 

► The maximum is taken at T = T^ax ~ 7c. We 
observe a sharp maximum of Ayi2 signalling that 
the asymmetry of the gluonic medium is peaked 
around the phase transition. 

► There is a symmetric point, T — Tq ~ ^Tc, where 
the asymmetry vanishes. This point is realized in 
the deconfinement phase sufficiently far from the 
phase transition at a temperature approximately 
equal twice the critical temperature. 

The mentioned points divide the phase diagram into 
three separate regions: 

► Region 1: The confinement phase, T < Tc- The 
asymmetry is a positive monotonically increasing 
function of temperature in the confinement region. 

► Region 2: The deconfinement phase at relatively 
low temperatures, Tc < T < Tq. Here the asym- 
metry is still positive valued and monotonically de- 
creasing function of temperature. The asymmetry 
vanishes at the temperature Tq, which we estimate 
below. 

► Region 3: The deconfinement phase at high tem- 
peratures, T > Tq. Here the asymmetry is negative 
valued and monotonically further decreasing as a 
function of temperature. 

As we speculate below in Section IIVI these regions 
are characterized each by a particular dynamics of the 
Abelian monopoles. These are singular configurations 
of the gluonic fields responsible for the confinement of 
color in the low-temperature phase. Region 1 corre- 
sponds to the phase, where the monopoles are condensed. 
At higher temperatures, in Region 2, the monopolc con- 
densate melts into a monopole liquid [29], whereas at 
even higher temperatures, in Region 3, the liquid of 
the Abelian monopoles is suggested to evaporate into a 
gaseous state [l^. The transition between Region 1 and 
Region 2 is a true phase transition (turning the monopole 
condensate into a monopole liquid), whereas the transi- 
tion between Region 2 and Region 3 (the evaporation of 
the monopole liquid) is suggested to be a broad crossover. 
The interaction with electric charges plays an important 
role in the dynamics of monopoles [30|, |31I | , such that we 
expect an effect of the changing monopole dynamics on 
the electric-magnetic asymmetry of the condensate. 

The deconfining transition at T = Tc has a noticeable 
effect on the asymmetry of the condensate since at this 
temperature the electric part of the condensate is max- 
imally dominating over the magnetic one. This domi- 
nance is rapidly decaying just above Tc- We define the 
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maximum as 



= mjix A^2 (T) = A^2 (Tn 



(35) 



We estimate the numerical value of the maximum as 

A^l''' = [2.26(4) TJ^ = [690(22) MeV]^ (36) 



and the temperature 



Tniax - 1.00(3) T, = 305(12) MeV. 



(37) 



The symmetric point is realized at a temperature at 
which the asymmetry of the condensate is vanishing. Our 
estimate of the symmetric point is 

To = 2.21(5)7; 675(23) MeV, Aa2(To) = 0. (38) 

The change of sign of A^2 happens in the dcconfinemcnt 
phase at a temperature approximately twice the decon- 
finement temperature. In order to accurately estimate 
the position of the symmetric point psp we performed a 
specific fit of the asymmetry throughout the deconfine- 
ment region. 

The high-temperature fit is done with the help of the 
following fitting function 



Afit 



.(0) 



[at high temperatures] (39) 



where A^2 and / are two fitting parameters. The form 
of the fitting function is inspired by the analytical ex- 
amples provided by the photodynamics (I19p and by the 
Abelian Higgs model (j27p in their high temperature lim- 
its. Surprisingly, the fit works very well not only at very 
high temperatures, but it also works down to the tem- 
peratures as low as 1.5 Tc. We obtain - with a Xdof ~ ^ 
- the following best fit parameters: 



^^2 — 



[0.894(14) T^f = [274(8) MeV]^ 



/ = 0.164(4), 



(40) 



The fit is shown in Figure [5] by the solid line. Accord- 
ing to the high-temperature limits in photodynamics (|19p 
and in the Abelian Higgs model ([?71) . one could have ex- 
pected that each color component of the gluon would 
give the same contribution 5^/12 to the coefficient / in 
Eq. (gni). For Ng = = 3 free gluons this coefficient 

should be equal to f — (?^/4, indicating that g ^ 1. This 
result is expected since in the considered temperature 
regime the theory is still in a strongly non-perturbative 
regime. One should note that the quadratic fit ([55)1 may 
mimic some other, nontrivial T dependence of the asym- 
metry, the exact form of which is difficult to figure out 
at our accuracy. The fit is convenient for the estimation 
of the symmetric point ([55)1 which follows from Eq. ([55)) 
and Eq. ([iil)). 

In order to emphasize the approach of the electric- 
magnetic asymmetry of the condensate to the asymptotic 
behavior at high temperatures we present in Figure [5] the 
asymmetry normalized by the temperature squared. 




T/T 

c 

FIG. 3: (Color online) The same as in Figure[2]but for electric- 
magnetic asymmetry of the j4^-condensate ([9]) normalized by 
the temperature T^. The inset zooms in the high-temperature 
region of the fit. The peaked value at T = Tc, Eq. ([40)) . and 
the asymptotic value at T ^ oo, Eq. (|36|l are indicated by 
vertical dotted and horizontal dash-dotted lines, respectively. 



The low-temperature limit is especially interesting in 
view of the existence of two different options: although 
Yang-Mills theory possesses the mass gap, the asymme- 
try would not necessarily be exponentially suppressed, 
as suggested in Eq. (jTUjl. Indeed, the asymmetry could 
behave polynomially, according to our calculation ()28p in 
the broken phase of the Abelian Higgs model. In order to 
figure out the behavior of the asymmetry we have made 
an interpolating fit, which includes both options: 



A fit 
A^2 



CaT^[y) exp{-m/T}, [at low T] . (41) 



The polynomial behavior would be realized if m — 0, 
while the exponential suppression is in effect with m ^ 0. 

We show the low-temperature asymmetry in Fig- 
ures Hl^a) and (b). In order to visualize the exponen- 
tial behavior in Figure [DJa) we took the logarithm of the 
asymmetry normalized by the dimensional factor T^ and 
then multiplied the logarithm by —T/Tc- If the behavior 
of the asymmetry is proportional to the Boltzmann-like 
exponential function (without the polynomial prefactor) 
then the data must be linear in the low-temperature re- 
gion. This is indeed the case in the region close to the 
phase transition while as temperature gets smaller then 
the deviation from the linear behavior becomes more vis- 
ible. Figure (D^b) corresponds to the different normaliza- 
tion factor under the logarithm, T^. 

First, we performed the fit using all three parameters 
Ca, V and m. The best fit parameters are presented in 
Table HI and the corresponding curve is shown in Fig- 
ure dj^a) by the dotted line. The fit was performed us- 
ing all available values of the asymmetry in the low- 
temperature region, T < Tc. As one can see from the 
Table, the exponent v is quite close to 2, so that we have 
fixed v = 2 and performed another fit. The quality of 
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both fits are characterized by almost the same value of 
X^/d.o.f., also presented in Table HI The fitting curve is 
shown in FigureHJa) by the solid line. Despite very good 
visual coincidence of the numerical data and the fitting 
curves, the high value of x^/d.o.f. « 6 is due to the fact 
that the data corresponding to different lattice volume 
are somewhat scattered with respect to each other in the 
region close to the phase transition. Since the transi- 
tion is of second order, we attribute the high value of the 
X^/d.o.f. parameter to finite- volume effects. 



V 




Ca 


xVd.o.f. 


m, MeV 


1/m, fm 


1.82(28) 


0.80(19) 


10.24(1.96) 


6.6 


244(58) 


0.80(19) 


2 [exact] 


0.66(2) 


9.00(32) 


6.2 


201(8) 


0.98(4) 


[exact] 


2.04(3) 


32.96(64) 


14. 


622(18) 


0.32(1) 


2.97(3) 


[exact] 


4.84(32) 


9.6 


[exact] 


oo 



TABLE L The best fit parameters v, m and Ca and the pa- 
rameter /A.o.l. describing the electric-magnetic asymmetry 
by the fitting function H4ip in the low-temperature region. 
The mentioned fitting curves are shown in Fig. [l] We also 
indicate the masses and the correlation lengths A = 1/m in 
physical units. 

The polynomial prefactor plays an important role. If 
we set 1/ = 0, the fitting function reduces to a purely 
exponential behavior. In this case the quality of the fit 
deteriorates drastically. The best fitting curve is shown 
in Figure HJa) by the dotted curve. Similarly, if we set 
the mass to zero, m = 0, and fit the data using only Ca 
and v as free parameters, the quality of the fit gets worse, 
represented by the dashed line in Figure [DJa). Moreover, 
according to Table U in this case the polynomial behavior 
would be - with a good accuracy - proportional to and 
not to the fourth power as one might have guessed from 
the example of the massive Abelian vector model ([28l) . 
Thus, the most plausible fit is characterized by an expo- 
nential suppression with a quadratic prefactor: 

^% = Ca exp{-TO/T} , (42) 

with 

Ca = 9.00(32), TO = 201 (8) MeV. (43) 

The best fit curve corresponding to the quadratic polyno- 
mial behavior ([1^ is also shown in FigurelDJb). We show 
the specific function of the data, —(T/Tc) log A^2/T^. 
This function must be linear if the data obey the law (|42l) . 
and this seems to fit the data. However, in order to figure 
out this fact with confidence one needs to calculate the 
electric-magnetic asymmetry at lower temperatures (the 
lowest temperature available to us in our simulation was 
T w OATc). 

There are two important remarks now in order. 

► Firstly, our data suggest that the behavior of the 
electric-magnetic asymmetry at low temperature is 
not proportional to the fourth power of tempera- 
ture, as one could guess from the simple case of a 



massive photon ((28)) . The behavior is rather ex- 
ponential (|32]) with a polynomial prefactor. Ac- 
cording to our discussion at the end of Section [TTbI 
a purely polynomial behavior is guaranteed if the 
propagator in momentum space tends to a cer- 
tain non-vanishing limit with vanishing momen- 
tum. The propagator in the Abelian Higgs model 
provides us with a clear example of such a behavior. 
According to our discussion, in the Abelian Higgs 
model the very reason of the purely polynomial 
behavior is the appearance of the "longitudinal" 
massless pole, in the integrand (jB2p of the in- 

tegral representation of the asymmetry (jAlip . The 
l/p^ term enters this representation in the com- 
bination D{p'^)/p'^, and if the propagator, D{p^), 
would vanish at low momenta as, for example, 
D{p^) oc p^, then the polynomial behavior of the 
asymmetry would change to an exponential one. 

Thus, the low-temperature behavior of the electric- 
magnetic asymmetry of the dimension-2 conden- 
sate in Yang-Mills theory may shed some light on 
the low-momentum behavior of the gluon propaga- 
tor in the momentum space. 

Indeed, the integral representation of the A^- 
condensate in the Abelian Higgs model (jAlip is 
similar to the one in Yang-Mills theory (|32p . such 
that the same considerations may apply. The expo- 
nential suppression of the condensate at low tem- 
peratures may signal that the thermal gluon propa- 
gator at low momenta behaves softer than the tree- 
level propagator. We should admit that we do not 
have data at low enough temperatures to prove this 
fact firmly. 

► The second interesting observation is that the mass 
which governs the exponential falloff of the asym- 
metry at low temperatures is not the glueball mass, 
as naively expected ([TT]). According to Eq. (^5]) 
this massive parameter is much smaller than the 
glueball mass (jS]), namely of the order of Aqcd, in 
other words, of the order of the critical tempera- 
ture, Tc, Eq. Thus, the characteristic length 
that describes non-perturbative effects related to 
the A^-condensates may be as large as one fcrmi. 
This could explain the fact that the exploration of 
the low-momentum asymptotics of the gluo n prop- 
agator requires relatively large lattices [32| . 



IV. DISCUSSION: CONFINEMENT AND 
ASYMMETRY OF CONDENSATE 

It is extremely interesting to understand a physical 
reasons behind the observed behavior of the electric- 
magnetic asymmetry. Figure [21 First of all it is strik- 
ing that in the confining region the condensate is not 
electric-magnetic symmetric. Instead, the asymmetry is 



10 




(a) (b) 

FIG. 4: (Color online) The behavior of the electric-magnetic asymmetry in the low-temperature region. The lines represent 
the fits by the function (|4ip which are discussed in the text. Notice the different normalization of the expression under the 
logarithm: in (a) the normalization factor is 1/T'c while in (b) this factor is l/T^. 



positive and is growing with temperature up to its max- 
imum which is realized just at the deconfinement phase 
transition. We know from the considered example of the 
Abelian Higgs model that a finite mass gap alone could 
not result in a positive value of the asymmetry. This sug- 
gests that the asymmetry of the dimcnsion-2 condensate 
ought to be related to the purely confining properties of 
the system. This observation is in agreement with the 
original suggestion made in Refs. 0, [s], as well as with 
numerical results of Ref. [ssj , in which a relation between 
the dimension-2 condensates and the confining string was 
discussed. Thus, one can conclude that color confine- 
ment (and its agents) may contribute to the unexpected 
behavior of asymmetry in the confinement phase. 

It is generally accepted that the confinement of color 
can be explained by the dynamics of either monopole- 
like gluonic configurations (for a review see Ref. (l3 |) 
or by string-like vortex configurations (a review can be 
found in Ref. [3]). In fact, these objects turn out to 
be interrelated physically [SJ, [H, [H, and geometri- 
cally ,1S] . Confinement of color in the low-temperature 
phase is caused by condensation of monopoles and by the 
spatial percolation of vortices. 

Since the asymmetry of the condensate is pre- 
sumably related to the color confinement one should be 
able to trace the asymmetry back to the dynamics of 
the monopoles (and of the center vortices, but below 
we discuss the magnetic monopoles only). Moreover, 
the interactions of the electric charges and the magnetic 
monopoles are imp ortant to determine the state of the 
monopoles [SO*, '315. Therefore, it is very natural to sug- 
gest that the monopole dynamics leaves its foot-prints 
in the electric-magnetic asymmetry of the dimension-2 
condensate. 

We suggest that the increase of the electric-magnetic 



asymmetry in the confinement phase is related to a grad- 
ual decrease of the monopole condensate, as witnessed 
also by the decreasing string tension with T — > Tc in 
SU{2) gluodynamics. Once the condensate has disap- 
peared at r = Tc, the asynimetry starts to become 
weaker. According to Refs. [2^[33,[3lj in this region the 
monopoles form a monopole liquid. The heating of the 
liquid leads to its complete evaporation only at higher 
temperatures. According to Ref. [1^ this should hap- 
pen around T « 2Tc which is pretty close to the point 
To where the asymmetry vanishes (j38l) . Thus, one can 
suggest that at T = Tq the liquid predominantly turns 
into a monopole gas. The negative value of the electric- 
magnetic asymmetry of the condensate is character- 
ized by the gaseous phase of the magnetic monopoles. 

V. CONCLUSION 

We have studied the electric-magnetic asymmetry of 
the condensate in the Landau gauge in three gauge 
theories at finite temperature in four space-time dimen- 
sions. 

► In photodynamics - the free theory of a massless 
Abelian gauge field - the asymmetry can be 
computed analytically. The asymmetry turns out 
to be negative for all temperatures being propor- 
tional to the temperature squared 

► In the Abelian Higgs theory (PO)) . for the London 
limit the analytical expression for the asymmetry 
is given by Eqs. ^ and jM]). This theory 
corresponds basically to a free theory of a massive 
gauge field. The asymmetry - shown in Figure [T]- 
is also negative for all temperatures, and the high- 
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and low-temperature limits can be computed, re- 
spectively, in Eqs. (^7)) and 



► In the SU (2) gauge theory we compute the electric- 
magnetic asymmetry using numerical simulations 
on the lattice. Figure [51 The temperature de- 
pendence of the asymmetry turns out to be unex- 
pected: at low temperatures the asymmetry turns 
out to be positive. It grows with increasing tem- 
perature, reaching a maximum around the critical 
temperature T « Tc. In the deconfinement phase 
the asymmetry drops rapidly with increasing tem- 
perature. At T = To = 2.21(5)rc = 675(23) MeV 
the asymmetry vanishes, and it becomes negative 
at higher temperature. 

In a spirit of Ref. [1^ we suggest that the low tem- 
perature asymptotics of the condensate is re- 
lated to the low-momentum behavior of the gluon 
propagator. Our data at relatively low tempera- 
tures, OATc ^ T < Tc, indicate that the asymmetry 
is suppressed exponentially (j42p providing strong 
arguments in favor of the infrared suppression of 
the gluon propagator at low temperatures. 

The mass, which governs the exponential low- 
temperature suppression of the asymmetry (|42p , is 
unexpectedly much smaller than the mass of the 
glueball. We found m = 201(8) MeV, somewhat 
smaller than the deconfinement temperature, Tc, 
of the pure Yang-Mills theory. The correspond- 
ing correlations length, A = 1/m is of the order 
of one fermi. As a by-product, this result suggests 
that relatively large (many fermi in one direction) 
lattice volumes are needed to pin-down the low- 
momentum asymptotics of the gluon propagator. 

The electric-magnetic asymmetry is most probably re- 
lated to the changing dynamics of the confining gluonic 
configurations, the magnetic monopoles. In fact, the 
monopole dynamics may be foot-printed in the electric- 
magnetic asymmetry of the condensate because the elec- 
tric degrees of freedom affect the properties of the mag- 
netic monopoles [13, HH and vice versa. More in detail, 
the regions of (1) positive and growing, (2) positive and 
decreasing, and (3) negative asymmetry coincide with 
the regions in which the monopoles should, according to 
classification of Ref. [1^, (1) be condensed, (2) form a 
liquid state and (3) form a predominantly gaseous state, 
respectively. 



APPENDIX A: THE ASYMMETRY IN 
PHOTODYNAMICS 

The asymmetry of the condensate in the momentum 
representation is 



(Al) 



n 



d^p 



{2nr 
d^p 
(2^ 



1544 (P,W„) - -^D,,{p,UJn) 



where the momentum-dependent structure function of 
the asymmetry is 

1 ^ 

G'a2(p,w„) = D44(p,t^«) - Ai(p,w„). (A2) 

In order to proceed further we use the following trick. 
It is well known [l^, [sl] that sums over Matsubara fre- 
quencies (|16p of the form (|A1[) can generally be converted 
into integrals, 

T FiPo = i^n) - 2^ Jdpo Re T(po) 

H — : / dpo friPo) ReT(po) , Po = «P4, (A3) 

TTl J 



where e +0, the real part of the function F is defined 
as follows, 



ReF(po) 



1 



T(po) + F(-po) 



and the temperature-dependent function 
/t(po) 



ePo/T _ 1 



(A4) 



(A5) 



is the Bose-Einstcin distribution of a bosonic particle 
with the energy uj — pq at temperature T. Equation (|A3p 
is valid for any analytical function F{pq) provided it does 
not possess poles at imaginary po-axis. 

It is easy to notice that the first term on the right 
hand side of Eq. (jA3[) is temperature-independent, and 
therefore it represents the zero-temperature part of the 
sum. Since the summation over Matsubara frequencies 
is usually supplemented with integration over the spatial 
momentum p, the zero-temperature part can generally 
be divergent in the ultraviolet region. The second term 
is the temperature-dependent correction which is usually 
finite because of the exponential suppression of the ul- 
traviolet modes with po ^ 0. 

The sum (jAip over the asymmetry structure function 



A^2(T) = A^.(r = 0) + <5A^.(r), 



with 



d^p 
(2^ 



(A6) 



(A7) 



12 



ioo+e 



dpo 
27r 



-ioo-\-€ 



and 



Ca(p, ipo) + Ca(p, -jpo) 



d3p 



2i7 (27r)3 



(A8) 



27r 



Ca(p, ipo) + Ca{p, -ipo) fripo) 



is expressed through the simple function 
One can explicitly check that 

A^2(r = o) -0, 



(A9) 



(AlO) 



as we have expected. This result is intuitively clear be- 
cause at zero temperature the difference between spa- 
tial and temporal directions disappears and therefore 
the space-time asymmetry of any quantity must be zero. 
Thus, only the temperature-dependent part of the 
condensates contributes to the asymmetry. 

At non-zero temperature we represent the asymmetry 
as the sum of the residues 



d-^p 



(All) 



The residue in the case of an nth order pole at z = a 
is defined as 



res /(a) 

such that 
for n ~ 1 



1 



lim 



(n — 1)! z^a dz 



(z-arfiz) ,(A12) 



res/(z) 



lim 

z — >a 



{z - a)f(z) 



(A13) 



for 71 = 2 res/(z) = lim — (z-aff{z) . (A14) 

z^a dz L J 

The quantity in the integrand, 

2 (p' + 3p2) /j,(po) 



2 {Rb CA{p,ipQ) f{po] 



3(po-|p|)2bo + |p|)2 



has a double pole in the complex pq plane on the positive 
real axis at pa = |p|. Using Eq. (|A14p and then Eq. (|Alip 
we get Eq. ([T9| for the case of the free massless photon. 



APPENDIX B: THE ASYMMETRY IN THE 
ABELIAN HIGGS MODEL 



In the Abelian Higgs model the instead of Eq. (jA9[) we 



get 



Ca{p,P4:) 



3p1 



3 (p^ + m?) 



(Bl) 



and, therefore. 



2 (Re CA{p,ipo) f{po)) (B2) 
_2 jp^ + ipDMPo) 



Thus there are two single (n = 1) poles in the complex 
Pq plane on the positive real axis at po = |p| and pq — 
i/pM-m^. Using Eq. (|A13p we obtain 



2 (ReC^(p,zpo) f{po) 

Vp^+m^ MVp^ + m^) - IpI /t(IpI)' 



3to^ 



1 MVp^ + m2) 
3 \/p'^ + m? 



(B3) 



Substituting this result in Eq. (jAlip we get Eqs. 
and 
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